Transfer and storage of qubits in the presence of decoherence 
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The effects of decoherence on the transfer and storage of coherent quantum states in hybrid 
systems are studied within the Caldeira-Leggett approach. In general, we find that a high transfer 
fidelity can be achieved even if the decoherence time is less than an order of magnitude larger than 
the transfer time, which is of order half a Rabi period and determined by the qubit-qubit coupling 
strength. We finally apply our results to assess the feasibility of a novel hybrid quantum memory 
system, comprised of the hyperfine qubit states of an ultracold atomic Bose-Einstein condensate 
and the flux qubits of a SQUID. 
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The implementation of a quantum analog of random- 
access memory for qubits, qRAM, is currently one of the 
major goals of quantum information processing [1]. As 
with their classical counterpart, which can be found in all 
modern-day computers, the qRAM should be both fast 
and robust. Ideally, arbitrary coherent quantum states 
can be quickly passed to, stored, and retrieved from the 
qRAM with high fidelity. Unlike in classical systems, 
these processes have to be performed with an unavoidable 
continuous loss of information, caused by the detrimen- 
tal influence of decoherence. The minimization of the 
effects of decoherence on transfer and storage of quan- 
tum information is therefore at the forefront of imple- 
menting a qRAM architecture and of quantum comput- 
ing in general. To realize this demanding goal, quantum- 
hybrid systems represent some of the most interesting 
and promising candidates [2-7]. 

Within a hybrid architecture one subsystem, or qubit, 
works as part of the quantum central processing unit 
(CPU), while the other component operates as the 
qRAM. One of the major benefits of combining seemingly 
disparate systems is that they can each reinforce various 
advantages, while compensating for deficiencies, of either 
subsystem. For example, recent experiments involving 
the nuclear and electron spin states of nitrogen- vacancy 
(NV) centers in diamond, and a sophisticated dynamical 
decoupling of qubits, have produced record decoherence 
times for a solid-state system [8], on the order of a sec- 
ond. Along with constructing a qRAM architecture from 
either a single NV center [9] or an ensemble of such states 
[10, 11], a variety of other hybrid systems have been ex- 
perimentally realized, such as superconducting qubits or 
atomic / molecular quantum gases coupled to electromag- 
netic and nanomechanical resonators [12-16]. 

By their very nature, each component of a hybrid 
system can have vastly different coherence times. Fre- 
quently, as is the case in many of the previous examples, 
the CPU element originates from a solid-state system, 
where the qubit states, |0) and |1), span the low-energy 
Hilbert space of some effective quasiparticle. While solid- 



state qubits have many advantages, such as scalability 
and controllability, they typically can have small coher- 
ence times. This is a major hurdle for quantum infor- 
mation processing. In contrast, the memory qubits may 
not be as easily constructed on a large scale, but they 
potentially have much longer decoherence times, making 
them ideal for storage of quantum information. 

The transfer of a quantum state from one qubit to the 
other is commonly done in experimentally realized hy- 
brids via the Rabi process, cf. e.g. [8, 13], dynamically 
controlling one qubit's level spacing, bringing the two 
subsystems into and out of resonance for half a Rabi cy- 
cle. The Rabi period is determined by the qubit-qubit in- 
teraction. Ideally this coupling should be strong enough 
to enable a fast transfer, but weak enough to not excite 
states outside of the qubit Hilbert space. To the best of 
our knowledge, a general study on the effects of decoher- 
ence during qRAM memory operations, i.e., state trans- 
fers, is however missing. In what follows, we undertake 
such a study, determining a practically applicable lower 
bound of the transfer fidelity which can be obtained as 
a function of the ratio of decoherence to transfer time. 
This bound is of fundamental importance for the devel- 
opment of current and future hybrids, when the latter 
ratio becomes of order unity. 

We simulate the transfer and storage of a state from 
one qubit to another in the presence of decoherence and 
determine the transfer fidelity. Within the Caldeira- 
Leggett model of dissipation and generalized spin-boson 
coupling [17] considered here, we find the limiting deco- 
herence time is determined by the energy-relaxation time 
T\ . Naturally, one would expect that to obtain a high fi- 
delity transfer, the transfer time should be much smaller 
than the coherence time. However, we find that to obtain 
a fidelity greater than 95%, a T\ time of only about three 
times larger than the qubit-qubit Rabi period, which in 
turn is larger than the transfer time by a factor between 
two and five (see also below) , is required. As a direct ap- 
plication of these results to a concrete system, we propose 
a novel hybrid qRAM system, composed of a flux qubit 
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(SQUID) and the hyperfinc states of a trapped atomic 
Bosc-Einstcin condensate (BEC). With current technol- 
ogy, the relevant time scales of decoherence and Rabi 
period can be of the same order. Therefore, the inclusion 
of decoherence effects in a state transfer simulation are 
needed to verify the feasibility of this qRAM architecture. 

For the coupled qubit-qubit system, the initial state is 
prepared in a parent qubit (p) and passed to a memory 
qubit (m). The Hamiltonian is taken to be (h = 1) 

H q - q (t) = — ^ e -^o-, - (o. )®<r x , (i) 

where A@B = A®t + t(£)B, Oi arc Pauli matrices, and 
the time dependence of the level spacing of the parent 
qubit, u) p (t) = uj m W(t). is used to bring the two qubits 
into and out of resonance [18]. The Rabi coupling f2 is 
assumed complex, due to the symmetrical a x ® <r x and 
Cy ® & x interaction terms. Although the form of the 
qubit-qubit coupling is chosen to coincide with the hybrid 
system proposed below, the following method and results 
can be easily generalized to include other forms. 

There can be many sources of decoherence, originat- 
ing both internally and externally. These sources will be 
collectively modeled phenomenologically, by the common 
Caldeira-Leggett type method of coupling the qubit's de- 
grees of freedom to a bosonic bath [17]. Assuming the 
decoherence of the parent qubit is the dominant limiting 
factor on state transfer, such that the decoherence of the 
memory qubit can be neglected in comparison, we only 
consider the parent qubit to be directly coupled a bath. 
The qubit-bath coupling is taken to be 

V = Y, X k^®(h + bl), (2) 

k 

where a = {<J Xl a y ,a z ) and A fc = (A^,A^,A|) is the cou- 
pling constant to each bath mode, labeled by k. The 
total Hamiltonian of the qubit-qubit-bath system is then 

H(t)=H cl . (l (t)®'£e k btb k + l®V, (3) 

k 

where determines the dispersion of the bath bosons. 

To accurately determine the fidelity of a state trans- 
ferred in the presence of decoherence, we numerically 
solve the time-dependent master equation for the reduced 
density matrix of the two-qubit system p{t). If at time 
t = the density matrix of the qubit-qubit-bath sys- 
tem is given by p(0) = p(0) <g> e _/ ^ b /Z h , with (3 = 
(feeT) -1 (the inverse temperature), Hh = Sfc e fc^fc^fcj 
and Z\j = Tre"* 3 ^, then the full density matrix at a 
later time is given by p(t) = U H (t)p(0)u}j(t), where the 
time-evolution operator Uait) is the time-order exponen- 
tial of the full Hamiltonian (3). Changing to the inter- 
action representation, p\(t) = W H (t)p(t)UH (t), where 
H Q (t) = H q . q (t)(BH h and U Ho (t) = Uh^)^"^ ■ 



Within the Markov approximation [19] the equation of 
motion can be written 

d t Mt) = -»[yi(t),ft(o)] - f dt' [y x (i), [ft &(*)]]■ 

Jo 

(4) 

In addition, we also make the Born approximation, which 
neglects correlations that develop between the qubit sub- 
system and the bath. In other words, at later times 
the density matrix is taken to be of the form p{t) ~ 
p{t) ® e~@ Hh /Zb. An equation of motion for the reduced 
density matrix can then be found from (4) by tracing 
over the bath; pi(t) = Ti'bPi(i). Assuming an isotropic 
qubit-bath coupling = Afc(e x + e y + e x ) = Afcn and 

Trbe~^^ b 6^ = 0, the equation of motion for the reduced 
density matrix is given by the equation 

d t p l (t)=- [ dt' [ deJ(e){coth(ej8/2)cos(e(*-t')) 
Jo Jo 

x [l®n-<ri(i), [l<g>n-<Ti(£'),/9i(0]] - % sin(e(i - t')) 
x [ign-ai^llgTO-oiCO,/?^*)}]}, (5) 

where the time dependence of the operators is now given 
by and J(e) = Efc A fc<5(e - e fe ) is a coupling 

weighted density of states or spectral density of the bath. 
A common choice of J(e) is the so-called ohmic bath, with 
J(e) oc e up to some upper cutoff u c . The form used 
here is J(e) = 7yee~ c / Wc , where rj is the dimensionless 
dissipation or friction constant. The energy integral over 
e can be done analytically, but is not important for the 
following results. 

Next, we numerically solve the 16 coupled differen- 
tial equations of (5) for the reduced density matrix el- 
ements using a fourth-order Rungc-Kutta method with 
a time step of to m t = 1 and parameters: u> c = 100 w m , 
j5 = lOOw^ 1 , |0| = 0.01 u m . Additionally, the nontriv- 
ial time-evolution operator Uh _ (t), which determines 
the time dependence of the operators in (5), and the fi- 
nal time integral have to be numerically found at each 
Runge-Kutta time step. The Rabi frequency used for 
the simulations allows for the calculations to be done in 
a reasonable time frame, while still within the weak cou- 
pling regime. Furthermore, the physically relevant quan- 
tity is the dimensionless ratio of the decoherence time to 
the time needed to transfer a state, approximately half a 
Rabi period Tr/ 2 = ■ 

Figure 1 shows the obtained decoherence times of the 
parent qubit as a function of the ohmic friction constant 
r\. These were found by first decoupling the parent qubit 
from the memory one, by setting Q = in (1). At 
time t = the system is then put in the pure state 
|*) = 2~ 1 / 2 (|00) + |01)), where the computational ba- 
sis states \ij) = \i) m (8 |j) p are the energy eigenstates 
of the uncoupled system. The energy-relaxation time 
T\ and phase-dccohcrcncc time T 2 arc then found by 
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FIG. 1. The calculated decoherence times of the parent qubit 
are shown in units of half the Rabi period T R / 2 = 2pJ' 
with |fi| = 0.01 com, as a function of the qubit-bath cou- 
pling. At time t = the system is placed into the state 
j*) = 2~ 1/2 (|00) + |01)) and then time evolved according to 
Eq. (5). The fit to the Ti data is given by f(x) = 5 - 10 * 10 ~ 4 

and for Tb by f(x) = The top panel shows an exam- 

ple of the decay of the elements of the density matrix as a 
function of time. Ti is found by fitting the decay to a func- 
tion of the form coe _t//Tl , while T 2 is extracted by fitting the 
off-diagonal decay to cq cos(cit)e _t ' /T2 . 



time dependent but decoupled qubit-qubit Hamiltonian 
H q . q in Eq. (1). 

Figure 2 shows the obtained fidelities as a function of 
Ti times, for a direct |01) — > 1 10) and superposition state 
transfer. As one can see, the fidelity is strongly depen- 
dent on the initial state. In particular to transfer the su- 
perposition state 2" 1 / 2 (|00)+i|01)) -> 2- 1 / 2 (|00)+i|10)) 
a fidelity greater than 90% can readily be achieved for 
Ti ~ T R/2 ; while for the transfer |01) [10) a fidelity 
of only about 50% is found. This is primarily due to 
the fact that for the superposition state the initial wave 
function has already a significant overlap with the final 
transferred one. Thus, even before any transfer steps are 
taken, it starts with a 50% fidelity. Generally, we find 
that to obtain a fidelity of at least 95%, Ti should be no 
smaller than about three times the Rabi period. Addi- 
tionally, because of the bath, after the transfer the state 
can still slowly dissipate from the memory qubit — even 
when the two subsystems are out of resonance. This oc- 
curs on a time scale on the order of (w m /|f2|) 2 Ti. This 
finite lifetime comes from a process where the memory 
qubit state |1) can relax, exciting the parent qubit's state 
from its ground state, which then decays into the bath. 
In diagrammatic many-body theory, this comes from a 
forth-order term in the expansion of the memory qubit's 



fitting the decay of the diagonal (Ti) and off-diagonal 
(T2) components of Pi(t) to a decaying exponential of 
the form e~*/ Ti [20]. As can be seen from Fig. 1, within 
this model T2 « 2Ti; therefore, the T\ energy relaxation 
time is the major limiting factor of any state transfer. 
The Ti time can also be found by calculating the life- 
time of the excited state r from diagrammatic many- 
body perturbation theory. The lifetime is then defined 
as t _1 = — 2Im£(u;p) [21], where £(w) is the parent 
qubit's self-energy. To second-order in the qubit-bath 
coupling, which contains the leading term describing the 
relaxation of the qubit and excitation of a boson in the 
bath, Ti ~ t = (4nriujp)~ 1 e u ' lp / ulc , which agrees well with 
the calculated values of the above model in Fig. 1. 

To transfer a quantum state from the parent qubit to 
the other, one first initializes the system by bringing the 
two qubits' level spacings out of resonance. Then, one 
prepares the parent qubit in an arbitrary initial state 
|*i) = a|00) + /?|01). Next, the two subsystems are 
brought into resonance for half a Rabi period T R y 2 - Fi~ 
nally, one takes the two qubits out of resonance, leaving 
the parent qubit in its ground state and its original state 
in the memory qubit |\&f) = ct|00) + /? 1 10) . The trans- 
fer fidelity is defined to be F(t) = ^(^ { (t)\p(t)\^{(t)), 
where p(t) is the reduced density matrix of the qubit 
system transformed back to the Heisenberg picture, and 
\^i(t)) is the final target state, time evolved under the 
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FIG. 2. (color online) The left column shows the fidelities for 
the direct transfer 1 1) — > |10), for various Ti times. Tr/2 is 
half a Rabi period, the approximate time required to transfer 
a state. The solid line in the top plots is the time profile of 
W(t) [18], the function used to bring the two qubits into and 
out of resonance, while the solid line on the bottom plots is 
a guide to the eye. To obtain a fidelity of more than 90% for 
this direct state transfer requires a Ti > 5T R/ / 2 . The right 
column shows the fidelities for the transfer of a superposition 
state 2- 1/2 (|00)+i|01)) -)• 2~ 1/2 (|00) +i|10)) , for different Ti 
times. Similar to the direct state transfer, to obtain a fidelity 
of more than 90% only requires a Ti > T R / 2 
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FIG. 3. A schematic representation of a SQUID-BEC hybrid. 
The flux qubit is created from superposing left- and right- 
going current states of the SQUID, coupled by its magnetic 
field to the hyperfine transition of an atomic BEC, where the 
Rabi process periodically flips one atom from a given hyper- 
fine spin state to the other. 

self-energy. Hence, we conclude that decoherence of the 
parent system has to be taken into account in general 
also after the qubit has been transferred to the qRAM. 

The immediate applicability of these results can be 
seen in the following novel quantum-hybrid qRAM ar- 
chitecture. The memory qubits consist of the hyperfine 
split ground state of a trapped ultracold atomic BEC, 
e.g, 87 Rb. These states have significant potential for a 
qRAM, as their T\ times are essentially infinite, on the 
order of seconds to minutes, longer than the lifetime of 
the BEC itself. Notwithstanding that currently there 
is no experimental data for their phase-coherence times, 
the use of this system for long-term storage of quantum 
information remains promising [22]. 

The BEC states are magnetically coupled to a flux 
qubit, which during a Rabi cycle excites a single atom 
from the macroscopically occupied ground state contain- 
ing N atoms (\) to an unoccupied excited state (f), see 
Fig. 3. The qubit subspace of the BEC is spanned by 

|0) B =(An)-V2(at)Vc), 
|l) B =[(iV-l)!]- 1 / 2 (4) Ar " 1 4|vac), (6) 

and the left- right current states of the SQUID [23], the 
uncoupled BEC-SQUID Hamiltonian can be written as 

Hbs = — o z © I ^ 2~~ °x ) , (?) 

where Whf is the hyperfine splitting of the BEC states (for 
87 Rb, Whf ~ 6.8 GHz), and e is the energy difference of 
the |L) and |R) current states with tunneling amplitude 
A(£), which can be dynamically modulated to bring the 
two systems into and out of resonance [24]. 

In the left-right current state basis of the SQUID, a 
current operator can be defined as follows; / ~ Ia z , 



where I = (L|i($)|L) ~ -(R|J($)|R) (exact equality 
holding when e = 0). The current operator is related 
to the flux operator by J = (<£ - & CX )/L with $ cx the 
classical externally applied flux and L the loop induc- 
tance. The magnetic vector potential of a loop of ra- 
dius R in the dipole approximation and polar coordinates 
reads A ~ ~n R 2 I sin6/(4r 2 ) a z [25]. The SQUID's 
magnetic field operator is then B = Vxi = B{r)a z . 

For an atom with total magnetic moment fi, (includ- 
ing orbital, nuclear, and electronic spin contributions) 
the magnetic dipole coupling of the BEC and SQUID 
is given by H- lnt = — J drrh(r) ■ B(r), where rh(r) is 
the second quantized magnetic moment density of the 
BEC; m(r) = *t(r)[/ik<r' **'(»•). In the qubit 

basis and neglecting Zeeman energy shifts, which sim- 
ply renormalize the energies of the flux qubit e — > e' , 
the coupling takes the form of that given in Eq. (1), with 
n = VNg f<l ■ /i tl4 ., where g a <7 , = J dr <j>* a {r)B{r)^ a , (r); 
here <f> a (r) is the center of mass atomic wave function for 
the trapped atoms. For a SQUID of radius 1 /im carry- 
ing a 1 mA current, and a center-to-ccnter BEC-SQUID 
separation of 50 /im, \g-\- 1 • Mt-ll ~ 100 H z - Therefore, 
given a typical BEC has TV = 10 6 atoms, Tr/2 is on the 
order of 10 //s, while current SQUIDs have decoherence 
times of the order of 1 fj,s [26] . 

Thus, from the results presented in the above, for the 
BEC to be a viable element of qRAM, the Rabi coupling, 
the SQUID decoherence rates, or some combination of 
both would have to be increased by a factor of approxi- 
mately 100, according to the above parameter estimates. 
The simplest way to achieve this would be to move the 
BEC closer to the SQUID. For instance, at a distance of 
10 yum the Rabi coupling would increase by almost two 
orders of magnitude. There are, hoewever, potentially 
detrimental effects caused by having the BEC in such 
close proximity to a (comparatively hot) surface, such as 
increased heating or spin flips caused by thermal emission 
from the SQUID [27, 28]. 

In summary, we find an effective quantum memory 
systems can be functional even if the time to trans- 
fer the state is of the same order as the smallest de- 
coherence time. We expect this to remain generally 
true, with only minor quantitative differences, for sys- 
tems with different qubit-qubit or qubit-bath couplings 
than those considered here (for example dissipation into 
a sub-ohmic bath) , because the ratio of transfer to deco- 
herence times should be relatively insensitive to the de- 
tails of the model. A detailed confirmation of this is left 
for future work. As an interesting and currently feasible 
application of this, a hybrid BEC-SQUID qRAM system 
was introduced. We have demonstrated that with the 
present coherence times of SQUIDs, the derived effects of 
decoherence on the transfer fidelity need to be taken into 
account to assess the functionality of the BEC-SQUID 
qRAM as a quantum memory device. 
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